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Abstract 

Everyday experience suggests that a ‘ruck’ forms when the two ends of a 
heavy carpet or rug are brought closer together. Classical analysis, how¬ 
ever, shows that the horizontal compressive force needed to create such a 
ruck should be inhnite. We show that this apparent paradox is due to the 
assumption of inextensibility of the rug. By accounting for a hnite extensi¬ 
bility, we show that rucks appear with a hnite, non-zero end-shortening and 
conhrm our theoretical results with simple experiments. Finally, we note that 
the appropriate measure of extensibility, the stretchability, is in this case not 
determined purely by geometry, but incorporates the mechanics of the sheet. 
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1. Introduction 

Localized bumps (or rucks) in a carpet are an everyday annoyance but 
have also been used as an analogy to understand a plethora of physical phe¬ 
nomena involving the sliding of two solid bodies. Examples include disloca¬ 
tions in plastic deformation [1], Schallamach waves in rubber friction [a Em, 
and even slip pulses in earthquakes |5]. Recently, a number of authors have 
focussed on the motion of such rucks, which can occur either rapidly (as 
when a rug is quickly shaken at one end) [HI Ej or slowly (as when a ruck 
‘falls’ down an inclined plane) [HI |9]. Surprisingly, however, some features 
of a static ruck remain poorly understood, including the conditions under 
which they form. 

At a superhcial level, the formation, or birth, of a ruck in a rug is similar 
to the buckling of the classic Euler elastica in- the two ends of a rug are 
brought closer together by a distance Af and a ruck forms to accommodate 
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the resulting excess length. Unlike the Euler elastica, however, the amount 
of buckled material is not equal to the whole system size — some material 
remains in contact with the ‘floor’ and the arc-length of the ruck, I, is not 
known a priori. 

Everyday experience suggests that I depends on the end-end compression 
Ai. To determine this relationship at the scaling level we follow Kolinski et 
al. |H]. For small rucks, simple geometry suggests that the height of the ruck 
d ~ For a rug of density ps, thickness t and bending rigidity S, 

we expect that the gravitational energy of the ruck ~ pspt x d x I, while 
the bending energy ~ Bid/t^Y x 1. Balancing these energies and eliminating 
d ~ (/A£)i/2 we End that 


Ai ~ 



( 1 ) 


As expected, the width of the ruck grows with increasing end-end compres¬ 
sion. 

How much compressive force is required for the onset of rucking? A 
simple calculation j6] reveals that the compressive force T required to form 
a ruck of size I satishes 


TB 


^ 80.76, 


( 2 ) 


with B the bending rigidity. Note that (|^ is precisely the classical result for 
the Euler buckling of a rod of length I with clamped ends m- Using ([^ to 
eliminate I from ([^ in favor of Ai reveals that the buckling load 


AfB 


( 3 ) 


which is divergent in the limit of very small end-shortenings Ai. To form a 
ruck from a flat rug we must pass through arbitrarily small end-shortenings 
and so we are led to the paradoxical result that to do so requires an inhnite 
compressive force! 

This divergence in the compressive force needed to form a ruck has been 
known for almost thirty years na. Previous authors have attempted to 
explain it as a result of a breakdown of the linearized beam theory used 
to obtain Q and @ da nans] while dS] conhrmed the validity of the 
linearized approach and showed decisively that a critical buckling load does 
not exist for a perfect, inhnite, continuous horizontal heavy elastica. Here, 
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we show using a combination of theory and experiment that the paradox is 
resolved by incorporating the hnite extensibility of the material. 

A key assumption in the preceding discussion is that the object, the 
rug in this case, has a hxed arc length, i.e. that it is inextensible. While 
this is true for objects with arbitrarily small thickness, any real object has 
a hnite, if small thickness, and so is, to a certain extent, extensible. The 
effect of hnite extensibility on the buckling of the classical Elastica has been 
studied by a number of authors [IZIIIHIE]. These analyses reveal that the 
crucial parameter governing the importance of extensibility is the ratio of 
the thickness, t, and total length, Ltot, of the beam or, equivalently, the von 
Karman number, 7 = . However, the ehect of 7 on the buckling 

threshold and the post-buckling behaviour is, in fact fairly small [H]. An 
unusual exception is the vibration frequencies of some modes of buckled 
beams, which may be diherent for perfectly inextensible beams (7 = 0 ) and 
asymptotically inextensible beams (7 —)■ 0 ): the limit is singular [ 18 ] . 

In this context, the classic prediction that the compressive force required 
to give birth to a ruck diverges clearly suggests the possibility that the hnite 
extensibility must ultimately play a role in regularizing the problem; we 
expect that inclusion of this effect will lead to a hnite compressive force and a 
hnite ruck size at the onset of buckling. Nevertheless, this suggestion appears 
not to have been made previously [T2| and so we study this problem here. 
Using a combination of theory and tabletop experiments, we show that indeed 
the hnite extensibility of any physical (rather than idealized) heavy elastica 
does give rise to a hnite ruck height at birth and, hence, a hnite buckling 
load. This is similar to the problem of the growth of a rod inside a curved 
cylinder where it was recently shown that hnite extensibility regularizes the 
otherwise inhnite pressure that would be applied on the cylinder prior to 
buckling I2D1. However, unlike this related problem, and previous studies of 
the role of extensibility, we hnd that the relevant parameter governing the 
role of extensibility is not simply determined by the geometrical properties 
of the object (namely t and Ttot) but depends also on its material properties 
(namely its Young’s modulus E and density Ps). 

2. Theoretical analysis 

2.1. Problem formulation 

We model the rug as a heavy elastic sheet of solid density ps, length 
thickness t and width b. This sheet rests on a solid horizontal surface, which 


3 



Figure 1: (a) A rug placed on a floor often ends up with a bump (a ruck) in it caused by the 
relative motion of its ends, (b) A laboratory ruck formed using sheets of Polyvinylsiloxane. 
For scale, the horizontal dimension of the square blue plate here is 15 cm. (c) Schematic 
sketch of the model problem considered here. A heavy sheet of natural length £oo, resting 
on a horizontal, frictionless surface, forms a ruck of lateral size I and height d when its are 
displaced by a relative distance A£. 
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we assume to be frictionless for simplicity, and is compressed by bringing 
the two ends a horizontal distance A£ closer together (see £gure[^). The 
sheet can respond to this compression in two ways: by compressing along its 
length and by buckling out of the plane. If the sheet buckles, it does so over 
a hnite region of width I < because the weight per unit area of the sheet 
{.PsQi) opposes the whole sheet losing contact with the surface. Assuming 
small transverse displacements, the prohle of the sheet y = w{x) satishes the 
linearized heavy elastica equation mm 


B—+T— 

dx'^ 


-pgt 


( 4 ) 


for |x| < 1/2. For 1/2 < |a;| < ioo/2 we have tc = 0: the sheet is in contact 
with the substrate. 

In Q, 5 = Ety[12{l - z/2)] is the bending stiffness of the sheet per 
unit width, with u the Poisson Ratio. As the problem is symmetric about 
X = 0, we consider only 0 < x < £ 00 / 2 , for simplicity. This gives rise to the 
symmetry boundary conditions 


n;'(0) = n;'"(0) = 0. 


( 5 ) 


Since w = 0for//2 < |x| < ioo/2, continuity of vertical displacement, 
together with force and torque balance ad x = 1/2 then give 

w{l/2) = w\l/2) = w"{l/2) = 0. (6) 

The system (|^-([^ is ostensibly a fourth order system with hve bound¬ 
ary conditions. However, the size of the buckled region, I, is unknown as is 
the imposed compressive force T (since we are assuming that the sheet is 
subject to a known displacement of its ends). In fact, then, we have a sixth 
order problem with the final condition coming from a relationship between 
compression and tension. Ordinarily, one might be tempted to specify that 
the imposed compression should be accounted for by the out of plane dis- 
placement, so that Ai = {w'Y dx. However, this condition implicitly 
assumes that the sheet is inextensible. Instead we return to Hooke’s law, 
which relates T to the strain, e^x within the sheet 

^ Et 

^ = “1-2®^ 

1 — 
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Since the strain can be expressed in terms of the horizontal and ont-of-plane 
displacements u and w, respectively, we then have 


Et 


du 1 / dtc 
dx 2 \dx 


Integrating from 0 to £ 00 / 2 , we hnd that 


^ 00(1 


froo/2 


2Et 


-T = - 
2 


dtc V 


A/ 




( 8 ) 


( 9 ) 


where we have used M(f'oo/2) = —^1/2 as the imposed displacement. (In 
the above derivation we have used that the compressive force T is constant; 
Hooke’s law ([^ then shows that the strain e^x is also constant.) 

The ordinary differential equation (|^ together with boundary conditions 
(§-(§) and the global constraint ([^ completely specify the problem. To 
make further progress, however, we hrst non-dimensionalize the problem us¬ 
ing the elasto-gravity length 

B 

- 10 

PsQtJ 

to rescale all lengths; ig emerges from the balance between bending and grav¬ 
ity that is expressed in (|^. We also introduce the dimensionless compressive 
force a = Tfg/B. 

The dimensionless equivalent of (|^ in the non-contact (buckled) region 

is 


IT"" + aW = -1. 

The boundary conditions (|^-(|^ are unchanged, while (|^ becomes 

r-Loc/2 /(jpp' 


S'crLoo = 


dX 


dX - AL, 


where {X, IT, L, Loo} = {x, tc, /, loo)/^g and 


12 £ 2 ’ 


( 11 ) 


( 12 ) 


(13) 


is the dimensionless ‘stretchability’ of the sheet. Note that while (13) has the 
form of the von-Karman number (7 = B/12E) discussed in the Introduction, 
the relevant horizontal length scale in the dehnition of S is neither the size of 
the sheet, nor the size of the ruck, I, but rather the elasto-gravity length. 


^g. 
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2.2. Analytical solution 

The governing equation (0 may be solved subject to the boundary condi¬ 
tions (|^-([^ to give the vertical displacement of the sheet in the non-contact 
region 


H'(A') = ^ 


cr^ 


cos(A/aX) 

cos{^/aL/2) 


+ 


{L/2)^-X^ 


Here, the dimensionless compressive force a must satisfy 

L^/a/2 = a 


(14) 


(16) 


where a ~ 4.49341 is the hrst root of tana: = x. Eqn (15) is the dimensionless 


version of ([^ while (14) is the classic shape of a horizontal heavy elastica, 
and has been derived previously PUS]. 

A relationship between the (unknown) compressive force, a, and the im¬ 


posed end-shortening, AL, is obtained by substituting (14) into (12) and 
yields 

= ^ + (16) 
Note that in the absence of extensibility. S' = 0, we recover the result P that 


L = f ^ 4 844 ^^i/7_ 


h 


( 17 ) 


which is the precise version of the scaling given in ([^. 

Equation (16) is the central result of our analysis. A sketch of AL{L) 
given by (16) with non-zero stretchability S' > 0 (see hgure [^) shows that 


for sufficiently large values of AL, two values of L may give a rucked solution; 
for small AL no such solution exists. Thus there is a minimum value of AL 
for which rucked states can exist. In particular, this turning point occurs at 


Lc = 




5-7 


1/9 


AL^^ = 


?16 


4 


1/9 


24.7V 


a^/\SL^y/\ 


(18) 

Rucked solutions exist only when the displacement AL > ALmin and, at this 
critical point, the ruck solution has width Lc. The critical buckling load is 


given by substituting this value of Lc into (15); crucially this load is hnite. 


The existence of an arbitrarily large buckling load for a ruck was based 
on the assumption that rucked solutions exist for arbitrarily small end- 
shortenings, as discussed around (|^. The existence of a minimum end¬ 
shortening ALmin for rucked solutions resolves this paradox but leaves the 
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question of precisely when a ruck should form unanswered: while a ruck is 
possible for AL > AL^in does a system necessarily attain that state as soon 
as AL > ALjnin? Furthermore, since two rucked solutions are possible for 
AL > ALmin, which of these will be observed in practice? 


2.3. Multiple states 

We have seen that for AL < AL^i^ no rucked solution exists and thus 
conclude that the the imposed deformation must be accommodated by com¬ 
pression. However, when AL > ALmin, there are three possible states that 
the sheet can adopt: two rucked states (with different values of L, see figure 
it) and one compressed. 

To determine whether spontaneous formation of a ruck is energetically 
favourable in comparison to the compressed state, we consider the mechanical 
energy of the system. The total energy of the system consists of bending, 
stretching and gravitational energy and may be written 


E 

~2 


re.oc/2 


[\Bwl 


+ 2^ 


-T) {u^ + \wl) + Psgtw] dx. 


(19) 


Here, the first term in the integral represents bending energy, the second the 
elastic strain energy and the third the gravitational potential energy. (In 
Appendix A we show that minimizing this energy subject to an imposed 


end-end displacement leads to the heavy elastica equation (|^.) 

Note that when there is no vertical deflection, the energy reduces to 


^flat _ lrp\p _ 1 Et Ai 

or, in dimensionless terms. 


( 20 ) 


t/a,, 1 AL" 

2B/4 *SL^' 


( 21 ) 


When considering the energy of the buckled state, it is important to 
note that the stretching energy is reduced (compared to the flat case) by a 
reduction in the compressive force. The contribution of the whole energy 
that comes from stretching is ff^°°^^(—T)exx/2 dx = (1 — iy^)T^£oo/(4Et). In 







dimensionless terms the energy of the buckled state becomes 


fh>uckled 

25/4 


or ^2 rL/2 

^^ + 1 C2W^xx + W)dX 

Aa^SL^ 7L^ 

I? ^ 384a2‘ 


( 22 ) 


The behaviour of the energy of the buckled solutions as the end-shortening 
increases is shown in hgure 0). We see that the buckled solution with the 
largest value of L has the lower energy of the two ruck solutions (as might be 
expected since it is the least curved). However, we also note that immediately 
after the buckled solution appears, it has a slightly higher energy than the 
unbuckled state and thus is only metastable. (We note that a cusped energy 
curve with a small portion of the buckled state that is metastable was also 
observed in a related study [20] •) Nevertheless, the length of this metastable 
branch decreases for decreasing SL^o, and as AL increases, the energy of 
this buckled state is considerably lower than that of the unbuckled, purely 
compressed state. Therefore, our model predicts that for stiff systems where 
the stretchability is a regularizing parameter {SL^o -C 1), the transition from 
flat to buckled solutions will occur close to ALmin and without hysteresis; this 
prediction will be validated in ^ 


2.4- Validity of approximations 

Our theoretical results have been predicated on the assumption that the 
deflection of the sheet are small and that, in particular, the slope of the sheet 
remains small. We have then shown that buckling only becomes possible for 
AL > ALmin ~ Now, the height of this buckled shape d ^ L‘^ 

AL^/^ while the width L AL^/^; hence the typical slope of the delaminated 

shape is d/L ~ AL^^"^. For this slope to be small at the onset of buckling 
(and hence for our approach to be valid) we therefore need that AL^JJ^ 1, 
i.e. S'Loo <C 1. We also need to ensure that the width of the ‘blister’ at onset 
is smaller than the total length of the sheet, which requires {SL^oY^^ 4 -^oo- 
However, we envisage that this will hold whenever S'Loo is sufficiently small 
and Loo > 1- 


3. Experiments 

To test the predictions of our model, we performed experiments with 
sheets with a range of thicknesses and Young’s moduli. The sheets were made 
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Figure 2: (a) The response diagram of the system as a function of the end-shortening, 
AL. We see that ruck solutions become possible only AL > AL^in{S). Results are shown 
for different values of the parameter SLao, as indicated, and the relationship for S' = 0, 
L « 4.844AL^/^ is also shown (black dotted curve), (b) A comparison of the energy of the 
buckled solutions (colored curves) and the purely compressed solution (solid black curve) 
as the end-shortening AL increases. We see that the buckled solution is not necessarily 
a global energy minimizer as soon as it appears, but becomes so very shortly after its 
appearance. Here results are shown for a variety of values of SLoo- Note that there are 
two buckled solutions possible for AL > ALj^in (see a); that with the larger value of L 
has the lower energy, and is the one that we expect to see experimentally. 
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E (kPa) 

t (mm) 

Ps (kg/m^) 

ig (cm) 


Symbol 

440 

4.8 

1220 

4.55 

9.3 X 10-^ 

A, V 

180 

5.0 

1220 

3.47 

1.7 X 10-3 

□ , ■, X 

94 

7.0 

1180 

3.51 

3.2 X 10-3 

o, • 


Table 1: Material properties for the three sets of experiments performed here, together 
with the symbol used to indicate that series of experiments in figure 


from Polyvinylsiloxane (PVS, supplied by Zhermack) of different stiffness; E 
was measured directly to be in the range 95 kPa < E < 440 kPa. The 
PVS sheets are made by casting the polymer melt in a rectangular tray and 
allowing it to cure. The width was not varied systematically but all results 
presented here had 10 cm < 6 < 14 cm while the length was varied in the 
range 23.3 cm < ioo ^ 29.3 cm. By varying the volume of polymer melt 
used, the thickness of the resulting sheet can easily be varied, but in each 
case is extremely uniform (to within 4%), as has been observed previously 
IZH; here we use sheets with thickness lying in the range 4.8 ± 0.2 mm < 
t < 7.7 ± 0.2 mm. Detailed values of the experimental parameters, and the 
values of S to which they correspond, are given in table 

The polymer sheets formed in this way can be made to be very soft and 
relatively thick. However, they are in general tacky, adhering to surfaces. 
They therefore do not slip easily past the substrate as has been assumed 
in our model. To overcome this problem, we adhered thin bands of acetate 
perpendicular to the direction of compression (as can be seen in hgure [^) 
allowing the sheet to slide smoothly over the substrate without changing its 
elastic properties noticeably. 

Figure shows the results of experimental measurements of the dimen¬ 
sionless ruck height D = d/ig as a function of the end-end compression AL. 
The ruck height is easier to measure experimentally than the region over 


which contact is lost, L; eqn (14) gives that 


D = 


32a4 


{a'^ + 2-2seca) ^ 0.00241^. 


These experiments were performed with three different values of the stretch- 
ability parameter S (as described in table and agree well with the predic¬ 
tions of the model presented earlier for the appropriate values of S. 

As expected on energetic grounds, the sheets jump to the ruck solution 
almost as soon as it becomes available and, further, it is the solution branch 
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0 0.2 0.4 0.6 0.8 1 

AL 

Figure 3: Experimental measurements of the dimensionless ruck height, D = d/ig, with the 
(dimensionless) end-end compression AL. Results are shown for increasing compression 
(closed symbols) and decreasing compression (open symbols) with a variety of different 
stretching stiffnesses: S = 0.00093 (triangles), S = 0.0017 (squares and stars) and S = 
0.0032 (circles). The corresponding theoretical predictions are shown by the solid curves. 
We also plot here the theoretical predictions in the limit of S' = 0: in this limit the linear 
theory predicts that D « 1.326AL^^^ (dashed black curve); the prediction from the 
geometrically nonlinear theory [6] is also shown (solid black curve). 
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with larger L that is selected. We also note that we do not observe hysteresis 
in this system: the results when the compression is gradually released (open 
symbols) collapse with those obtained as the compression is increased (closed 
symbols), to within experimental tolerances. This distinguishes this problem 
from apparently related problems such as the delamination of a stiff sheet 
from a soft foundation [22], where hysteresis is large enough to be easily 
observed. 

4. Conclusion 

We have studied the onset, or birth, of a ruck in a rug and shown that the 
hnite extensibility of the sheet means that, in reality, rucking can occur at 
hnite end-end compression, and with hnite compressive load. This resolves 
an issue with the standard, linearized but inextensible theory that has caused 
some confusion in the literature. 

Our analysis shows that rather than taking any given size, blisters have 
a preferred size, below which it is energetically preferrable to accommodate 
the imposed end-end displacement by compressing the material uniformly. 
In the delamination of a stiff sheet from a soft foundation, the existence of 
such a preferred size leads to the formation of a series of blisters as the end- 
end compression increases [ 22 ] • One might then naturally ask why does only 
one ruck form, rather than a whole series, as Ai increases beyond the onset of 
rucking? To answer this, we note that the energy of a ruck B(P/l^ ~ At/'^ 
— dividing a given end-end compression between N > 1 rucks leads to a 
total energy N X (A£/Ar)5/7 _ JV2/7, which is energetically unfavorable 
compared with a single ruck. 

We note that in other problems where non-zero extensibility has been 
found to play a role, the relevant stretchability parameter has been deter¬ 
mined purely from the geometry of the system. For example, the growth of 
a rod conhned within a circular tube is controlled by the ratio of the rod 
thickness to the cylinder radius 1201 , while the buckling, snap-through and 
‘ringing’ oscillation of beams and arches depends on the ratio of thickness to 
beam length [TRIITU] . In contrast, for the birth of a ruck in a rug, the relevant 
stretchability is S' = Rather than being purely geometrical, this S 

contains the mechanics of the problem via the elasto-gravity length scale £g. 
This is a consequence of the importance of two force scales in the problem 
(gravity and elasticity) and, while not generic, such balances may give rise to 
similar behaviour when extensibility is included in related systems [ 231123 ] • 
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Finally, we note that there are likely other ways in which the divergence 
of compressive stress might be regularized in this, and other, problems. For 
example, we expect that a hnite substrate compliance might lead to wrinkling 
of the combined system [25l |26] at small end-end compressions but become 
energetically unfavourable compared to rucking at larger compressions. 
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Appendix A. Energy and a variational principle 


Considering only half of the system, i.e. 0 < a; < f'oo/2, the functional 
describing the system consists of the bending, stretching and gravitational 
energies of the beam, together with a constraint associated with the end-end 
compression. This constraint is enforced by a Lagrange multiplier F so that 
the appropriate functional is 


r-^oo/2 


u = 


[\Bwl^ + \{-T){u^ + \wl) + psQtw] dx - F ^+m(£oo/2) 


(A.l) 

Here, the hrst term in the integral represents bending energy, the second the 
elastic strain energy and the third the gravitational potential energy. 

On hrst sight, it appears that the factor \Tw‘^ within the integral means 
that the coefficient of Twxx in the corresponding Euler-Lagrange equations 
will be 1/2, rather than unity in (|^. However, a more careful consideration 
reveals that the heavy elastica equation (|^ is, in fact, recovered. To see this, 
we hrst eliminate T from (A.l), using Hooke’s law ([^, and write M(f'oo/2) = 
Ux dx. Performing the variation with respect to u we hnd that 


F = 


Ft 


z/" 


{Ua 




)=-T. 


(A.2) 


At the same time, the variation with respect to w gives 


BWxxxx “t“ TWxx psQt 0, 


(A.3) 
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which is the expected heavy elastica equation, eqn Q. 

From this, we see that (19) gives the energy of the system and is consis¬ 
tent with the natural variational principle used to derive the heavy elastica 
equation (Q. 
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